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Abstract 

Recent work investigating the interaction of magnetic islands with micro-turbulence has uncov- 
ered the striking observation of large scale vortex modes forming within the island structure [W.A. 
Hornsby et al, Phys. Plasmas 17 092301 (2010)]. These electrostatic vortices are found to be the 
size of the island and are oscillatory. It is this oscillatory behaviour and the presence of turbulence 
that leads us to believe that the dynamics are related to the Geodesic Acoustic Mode (GAM), and 
it is this link that is investigated in this paper. 

Here we derive an equation for the GAM in the MHD limit, in the presence of a magnetic 
island modified three-dimensional axisymmetric geometry. The eigenvalues and eigenfunctions 
are calculated numerically and then utilised to analyse the dynamics of oscillatory large-scale 
electrostatic potential structures seen in both linear and non-linear gyro-kinetic simulations. 
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Magnetic islands, generated by the tearing mode, have the effect of breaking the ax- 
isymmetric properties of the equihbrium, and can have a detrimental effect on the plasma 
confinement due to the radial component of the magnetic field that they introduce 1-^. 

Recent work on the interaction of turbulence with large magnetic islands has uncovered 
the striking observation of large scale electrostatic vortices forming within the island sep- 
aratrix (See Fig. [2]). These meso-scale potential structures generate E x B flows around 



the island, similar to the zonal flows which act as regulators of turbulence Q-llO(, having 
the effect of tearing up radially extended electrostatic eddies, and thus reducing the radial 
transport of particles and heat. On the contrary, it is seen that these vortices can enhance 
the heat flow within the separatrix by up to 50% by acting as a convective cell, having a 
detrimental impact on heat confinement in a toroidal plasma lll-ll3| , but also have a sig- 
nificant effect on the radial pressure profile, which in turn can effect the bootstrap current 
profile which determines the stability of the Neoclassical Tearing mode (NTM) 14j-|l7| . 

Zonal flows in fusion devices are intrinsically linked to an oscillatory mode known as the 
Geodesic Acoustic Mode (GAM) 18|-|20[. These can be excited non-linearly by primary in- 



stabilities such as drift waves or by interaction with energetic particles in fusion plasmas |2l[|. 
GAMs are generated on closed flux surfaces when perturbations in the E x B flow couple to 
axisymmetric pressure perturbations, by way of the curvature in the magnetic field causing 
a compression, to produce oscillatory electrostatic modes. Geodesic Acoustic modes have 
been extensively observed in a variety of laboratory plasmas 22N26| . The GAM was initially 
identified in toroidal symmetric systems, however, recently both theory and experiment 
has extended to include helical systems 271. |28[|. It is thought that plasma compressibility, 
which is the cause of the GAM, can have a significant effect on the growth rate of tearing 
instabilities 29l |. 

It is the oscillatory electrostatic structure of the GAM, and the regular oscillation period 
of the vortex seen in nonlinear simulations (See black trace in Fig. [Hand 2D slices in Fig. |2]), 
that leads us to believe that similar physics is responsible for the oscillatory vortex structures 
seen inside magnetic islands and it is this observation that is the basis of this paper. 

The period of the GAM frequency is determined by the compressibility of the plasma 
and is closely related to the sound speed and to the shape of the flux surface. The radial 
component of the magnetic field introduced by the tearing mode, produces magnetic islands 
which are seperate confinement regions within the plasma. Within these new confinement 




FIG. 1. Time trace of the amplitude of the smallest radial wave- vector (sideband) as a function 
of time for an island with poloidal wave-vector, kepi = 0.05. There are two frequencies present, 
the fast, quickly damped Geodesic Acoustic Mode and the slower, island Geodesic mode, which is 
very slowly damped and of larger amplitude. The amplitudes of the modes have been adjusted so 
that a direct comparison can be made. Without the presence of turbulence the mode is initialised. 
When turbulence is present, the vortex mode is generated. 



regions the GAM period is likely to be highly modified due to their helical structure. 

The paper is structured as follows. In sections [T] and [Tll the model is outlined and the 
Eigenfunctions and frequency of the Geodesic Acoustic mode in the presence of a magnetic 
island structure are calculated and analysed in section UTTl 

Sections IIVI and |V] will outline Gyro-kinetic simulations to study these dynamic struc- 
tures and then a comparison is made with oscillations seen in fully nonlinear gyro-kinetic 
turbulence simulations with magnetic islands. 



I. MATHEMATICAL MODEL 

We begin by calculating the eigenvalue equation for the geodesic acoustic mode in the 
ideal MHD limit, and modify the theory to take into account the change in the flux surfaces 
due to the presence of a magnetic island. We assume here that the evolution of the island 
is significantly slower than oscillation time of the mode. The islands that we consider are 
large (an island half- width of w/pi = 24), and are treated as a static structure, rather than 



a dynamic mode and behave as a separate confinement region within the plasma 



30|. 




FIG. 2. (Color online) Normalized electrostatic potential {(pN = ^(l^/Tp^, w/ Pi = 24) in the plane 
perpendicular to the magnetic field (outboard mid-plane). Black lines represent the perturbed flux 
surfaces calculated from the total parallel vector potential. The presence of the island embedded in 
the turbulence not only generates flows around the island structure but also large scale electrostatic 
potential structures within the island which fluctuate in amplitude and sign. The top panel shows 
a vortex with a positive sign, while the lower panel shows that the vortex has flipped sign at a 
point later in the simulation. 

For the equilibrium magnetic field we utilise the axisymmetric, large aspect ratio toroidal 
geometry, with a further assumption that the flux surfaces are circular. Here, e = r/R is 
the inverse aspect ratio, where r is a minor and R is the major radius. 

The helical angle, associated with the magnetic island mode is defined as: 

^ = m9 — n( — ut (1) 

where 6 and ( are the poloidal and toroidal angles respectively, and m and n are the poloidal 
and toroidal mode-numbers, u is the island rotation frequency, which in this analysis is set 
to zero. 

The perturbation due to the magnetic island consists of a helical flux component, 

ip = ip cos ^ 

ip is treated as a constant in accordance with the well utilised constant 
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(2) 

—^jJ approximation 



With the island present, it is possible to construct modified fiux-surfaces with Q as the 



island fiux-surface label. This has the form 
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\L = — - = 2 cost (3 

Where tq is the minor radius at the rational surface on which the island is sited and w is 
the island half-width. 

We are interested, primarily, in modes with a small poloidal wave-vector in the electro- 
static limit, so we keep our analysis to the ideal MHD equations. 

Here we present the linearised MHD equations and closely follow the procedure to calcu- 
late the Geodesic acoustic mode dispersion relation that is originally outlined by Windsor 



et. al 



18| . The equations are as follows: 



i9v 

p— = J X B - 

ot 

dp „ - 
g^ + pV.v = 

V0 = V X B 

V- J = 

p-'% - + v.Vbp-^) = (4) 

Where the equations are the linearised momentum, continuity. Ohms law, quasi-neutrality 
condition and the equation of state, respectively. The tilde denotes a perturbed quantity, 
p, J and p denoting the mass density, current density and pressure respectively. B, and 
V are the magnetic fields, electrostatic potential and the plasma velocity. 7 = 5/3 denotes 
the adiabatic constant. 

The fiux surface label, fi, satisfies the magnetic differential equation, 

B ■ Vfi = 0, (5) 

when a magnetic island is present. With this in mind, the fluid velocity has three components 
and can be written as: 

^ Vfi . B X Vfi . B \ 
v = exp(-za;t) + + J (6) 

The flrst term being the velocity across the flux surfaces, second, the velocity in the binormal 
coordinate and the last term being the velocity directed along the magnetic fleld. From 



Ohm's law in Eq.( H]), it follows that the electrostatic potential is just a function of the flux 
surface label, fl. 

It can be shown that the term across the fiux-surfaces, is zero and = ^ which is a flux 
surface quantity. This reduces the equation of state to p = yp, the equilibrium pressure 
given by, p = pkBT/nii. 

Taking the flux surface average of the linearised momentum equation yields in the radial 
direction, 

= J / J "b^^"^*^' 

where J is the Jacobian, whereas the parallel component can be used to obtain, 

Vs = -i^B ■ Vp. (9) 

Substitution into the continuity equation we obtain the following Eigenvalue equation for 
oj"^, the squared mode frequency. 



tu'^p 
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JdS 



p — Jab 



B 



(10) 



The first term representing the effect of E-cross-B flows associated with the compres- 
sion caused by the geodesic curvature within the magnetic island, while the second term 
represents the motion of sound waves parallel to the magnetic field lines. 

Multiplying the continuity Eq.(|4]) with the complex conjugate of the mass density (p*) 
and taking the flux surface average, we arrive at an integral equation for the dispersion 



relation analogous to the one given in 



IP 



JdS 



j \pYJdS - 

p — JdS 



+ / \V\\p? JdS 



(11) 
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In the above derivation, finite gyro-radius effects have been neglected and as such does 
not provide any information about the structure of the vortex mode perpendicular to the 
perturbed flux-surfaces Q. This would require a kinetic or higher order approach 36|] and is 
beyond the scope of this paper. 



II. EIGENVALUE CALCULATION 

We calculate the Eigenvalues and Eigenf unctions numerically by writing Eq.( [TU]) in the 
form: 



2~ / IP 

W p = - — 



B X v^] • 



|vn|- 



JdS 



C 



where: 



C 



P JdS 



(12) 



(13) 



which can be written in the form of a generalised Eigenvalue equation and then solved 
using standard methods. 



In magnetic island geometry |32|-|35| . V|| is defined as: 



1 d 
Vii = 

" Rq de 



(14) 



and remembering that k\\ = —k 



Ls is defined as, Lg 



Rq/s 
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ijL~'^o) / Ls and also ke = nq/r = m/r and the shear length. 



35| 



fell 



wms 



T- 



V2qeR 



(15) 



Where vq is the radial coordinate of the rational surface of consideration, the negative sign 
is chosen when (r — tq) is positive and vice-versa. 

We note here that the flux surface integral is defined by firstly taking the integral over 



the poloidal angle 6 then by an integral over the 
A, this flux surface average can be written as|17l. 



ne 



ical angle ^. For an arbitrary function. 



(A) 



y4ro(l + ecos^ 

^/Q + COS^ 
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(16) 



Consider the term J JdS. Firstly we note that, we use the simple circular cross- 
section, axisymmetric model for tokamak equilibrium, with the approximations, Vr = 
r,V^ ~ f and VC ~ I- 

If we utilise: 



Vn = ^^'-;'h + !^sin^9-^sin^C (17) 
r R 



Performing the 6 integral, which removes the terms of order e. We assume a small aspect 
ratio (r < R) and also small island width in relation to the minor radius(ti' < r) and as such 
the integral reduces to: 



Vfi|\,„ 327r /(fi + cosO 



JdS = / \ ^' d^ (18) 

52 B^w^ J V^ + cose ^ ^ 

Finally we treat the term J p BxV^^.VB^ JdS, for brevity the full details of this calculation 

1pm 

"r7 



can be found in Appendix A. Utilising B = BtC + Bg9 + ^ sin^r, neglecting the effect of 



the island on the field strength, utilising: 

Bt = Bo/ (l + e cos 9) (19) 
= ^/(1 + ecos^) (20) 

and the major radius varying according to: 

R = Ro{l + ecose), (21) 

the compression term Bxv^^-vn -^^ written as: 



B X VB^ ■ Vn -2 /V^2 + cos^ 

W' = 1^1 ,,V-2 

/"m en \ ^ \ , , 

This equation is substituted into Eq. (fT2|) . which forms a generalised eigenvalue equation. 
Both positive and negative r — tq sides of the magnetic island are considered with periodic 
boundary conditions and the eigenfunctions and eigenvalues calculated. The results are 
discussed in next section. 

Far away from a magnetic island, where the perturbation of the magnetic flux surfaces is 
smallest, it is expected that the dispersion relation will return to the form of the standard 
Geodesic acoustic mode. 
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Taking the limit of large fl and large aspect ratio, e < 1. Taylor expanding then per- 
forming the integrals we obtain: 



-^-#(2 + ^). (23) 

From Fig. [H we see that the high frequency GAM is present in both the simulations with 
and without electrostatic turbulence, having an identical frequency in both. 



III. EIGENFUNCTION ANALYSIS 

Comparison of the island modified GAM dispersion with the result for toroidal circular 
flux surfaces, shows that the oscillations are significantly modified by the presence of a 
magnetic island. The oscillation has a longer oscillation period that the standard GAM by 
a factor of approximately (kgPi)^, which, is of the order of 10^^. Fig. [3] plots the density 
eigenfunction in the helical island direction, kg being the islands' poloidal wavevector, which 
appears whenever a derivative, Vfi, is performed. Plotted are the four lowest harmonics, 
higher harmonics exist but are neglected here. The inlay shows the function in the poloidal 
angle, which is sinusoidal in the same way as the standard GAM, but of significantly smaller 
amplitude. The function being dominated by the helical direction. 

The compression in this case is supplied by the variation of the magnetic field, B, as 
we travel around the magnetic island. A net compression exists when the Eigenfunction is 
a symmetric reflection between the outer and inner half of the island solutions, which can 
also be thought of as a symmetric solution in the radial coordinate (See Fig. |3]for a simple 
depiction). When a radially asymmetric solution exists the net compressive effect is zero 
and the solution represents a pure sound wave (e.g. Right hand panel of Fig. H]). 

Unlike the case of the normal GAM, when the inverse aspect ratio is set to zero we obtain 
solutions that are pure sound waves and the eigenvalue of the symmetric and anti-symmetric 
solutions are identical. With a finite aspect ratio these two value diverge as the compressive 
part becomes larger. 



The damping rate is determined by kinetic effects as showed by Hinton and Rosenbluth 371 . 



however we see from the time traces of the potential amplitude seen in Fig. [T] that the 



damping rate is significantly slower than the rate for the GAM. 
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FIG. 3. The non-trivial density Eigenfunctions in the direction of the helical angle, ^ at 17 = 0.95. 
Shown in inlay is the Eigenfunction in the poloidal angle. Plotted here is one half of the island 
eigenfunction (i.e. positive or negative radial direction). The curves are numbered according to 
increasing Eigenvalue. 



X-point. 




FIG. 4. Cartoon representing two configurations of density perturbation around magnetic flux 
surfaces in a magnetic island, (left) reflected, symmetric solution between inner and outer island 
regions, producing an up-down density asymmetry and a degree of compression, (right) zero com- 
pression perturbation where the outer solution is simply a copy of the inner solution, producing a 
radially asymmetric density perturbation. 



The collisionless damping rate of the GAM, normahsed to the ion transit frequency, has 



been shown to have the following form 



39|, 
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where u is the oscillation frequency and Ut is the ion transit frequency, estimated by Ut ~ 
k\\Vthi ~ kgwsvthi/q in the presence of an island. With some algebra it can be shown that 
the ratio of the damping rate to the transit frequency is, ^ ~ 0.2. Where the subscript I 
denotes the island modified GAM. 

With the same analysis; taking the ion transit time for the normal Geodesic Acoustic 
Mode as, k\\ = l/{Roq), we can calculate the damping rate as 7g ~ 0.36vthi/Ro- 

Finally, the ratio of the normal GAM damping rate to the island modified GAM damping 
rate is ^ ~ 15. . From this rough calculation we see that the damping rate of the island 
modified mode is significantly smaller than the damping rate of the GAM. 

Both modes are evident in the oscillations seen in the trace, the faster GAM, and the 
slower island oscillation. We note here that the faster oscillation is persistent in the turbu- 
lence simulation as it is being continuously excited by the turbulence, while in the 'linear' 
simulation only an initial perturbation is possible which damps accordingly. 



IV. GYRO-KINETIC FRAMEWORK 



The oscillatory vortex mode was observed in gyro-kinetic simulations studying the effect 
of magnetic islands on drift-wave turbulence. Since the GAM is a linear mode it should be 
possible to excite them without turbulence present. 

Here the GAM oscillation is studied using the gyro-kinetic framework, with numerical 
solutions obtained using a modified version of the gyro-kinetic fiux-tube code GKW of which 

n 

more details can be found in |40|. The delta-/ approximation is used, in which the distribu- 
tion function is split into a background F and a perturbed distribution /. The final equation 
for the perturbed distribution function /, for each species can be written in the form 

^ + (.,b + v,).V/ + v,.Vj-^— g^ = S. (25) 

where 5* is the source term which is determined by the background distribution function, 
fi is the magnetic moment, v\\ is the velocity along the magnetic field, B is the magnetic 
field strength, m and Z are the particle mass and charge number respectively. Here, g = 
f + {Ze/T)vii{A\\)FM is used to absorb the time derivative of the parallel vector potential 
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dA\\/dt which enters the equations through Amperes law. The background is assumed to be 
a Maxwelhan {Fm), with particle density {n) and temperature (T) 



n 



exp 



-'th 



(26) 



which determines the source term, neglecting temperature and density gradients becomes: 



S 



Z 

~[v\\h + WD\-y{(p)FM. 



(27) 



The thermal velocity fth = \/'^T/m, and the major radius (i?) are use to normalise the length 
and time scales. Using standard gyro-kinetic ordering, the length scale of perturbations 
along the field line (-RVy ~ 1) are significantly longer than those perpendicular to the field 
(-RV_L ~ 1/P*)- Here, p^, = Pi/R is the normalised ion Larmor radius (where pi = rriiVth/eB 
and Vth = \/W^Jrni). 

The velocities in Eq. f l25|) are from left to right: the parallel motion along the unperturbed 
field (f||b), the drift motion due to the inhomogeneous field (v/j), and the motion due to the 
perturbed electromagnetic field (v,^). The drift due to the inhomogeneous magnetic field 
can be written in the form 40| . 



1 



mvu 



B 



+ /U 



B X V5 

52 ' 



whereas the motion due to the perturbed electromagnetic field 

b X Vx 



(28) 



(29) 



is the combination of the E x B velocity (v^; = b x V{(f))/B) and the parallel motion along 
the perturbed field line (v^^ = — b x Vf y (Ay)/!?). These two effects are combined into one 
velocity through the definition of a new field x = (0) ^ "^11 (^||)- Here, the angled brackets 
denote gyro-averaged quantities. 

The electrostatic potential is calculated from the gyro-kinetic Poisson equation which in 
Fourier space is 



ZspURsp 2ttB / dv\\dpJo{k±psp)gsp + 

sp 



■'sp 



Trsp 

^2^2 _ 1 fcjm^g 



[r(&sp) - 1]0 



(30) 



where b = \miiTii{k±p^R^c{/ ZB'^Y = kx_ being the perpendicular wave-number 

and Jn are zeroth order Bessel functions of the first kind. 
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FIG. 5. (top) The square of the oscillation frequency for individual code runs where the value of 
kg Pi, and hence the toroidal size of the island, is varied, plotted against the square of the wave- 
vector. The dashed line represents a linear fit. This shows an exactly linear dependence. (Bottom) 
The calculated ratio of the Eigenvalue and (kgpi)'^ as a function of the perturbed flux label 1^ from 
the 0-point = — 1 to the separatrix, Q = 1. The agreement between the calculated gradient from 
the (top) figure, uj]^ / [kg pi)"^ = 3.7 (represented by the dashed line, here ujn = ^iVthi/R) compares 
well with the value at approximately the island separatrix, ~ 0.9. The code was run with 50 
points in both the ^ and 6 directions. 



GKW uses straight field line Hamada 4l| coordinates {s,C,,il)) where s is the coordinate 
along the magnetic field and C, is the generalised toroidal angle. For circular concentric 
surfaces, the transformation of poloidal and toroidal angle to these coordinates is given by 



40| (s, C) = (^/27r, [qO — 0]/27r). Assuming the winding of the magnetic field is resonant 
[q = m/n) in the centre of the computational domain [ipQ = vq/ Rq = e, where r is the radius 
of the magnetic surface, and Rq is the distance of the centre of the surface to the axis of 
symmetry) and expanding q up to first order in Aip {Aip being the radial distance from the 
resonant surface and here m is the poloidal mode number), q = m/n + Aipldq/dtp) then 
yields 

^11 = i|| exp[27rin(C - sdq/d^Aip)]. (31) 
13 



The wave vector of the island is k^pi = 27rnp^,. GKW uses a Fourier representation in the 
plane perpendicular to the magnetic field. The periodicity constraint on the torus shaped 
magnetic surface then dictates a relation between the radial and toroidal modes. 
The half width of the island is defined by 



w = 2Jqij/sRB, 



(32) 



(where s = {l/q)dq/dilj is the magnetic shear) and the perturbed magnetic fiux, is related 
to the perturbation of the parallel vector potential by the relation, 



-RAii 



(33) 



Full details of the numerical implementation of t 
interested reader can find them in these papers [isl, 



le magnetic island is omitted here, the 
3. 



V. RESULTS AND COMPARISON 

Presented here are the results from simulations, which keep the kinetic electron effects 
with the true mass ratio of a Deuterium plasma. While we are studying the linear response 
to a perturbation, due to the set-up of GKW, the code must be run non-linearly for the 
plasma to feel the effects of the modified field lines due to the magnetic island. This is 
because, due to numerical reasons, the parallel vector potential of the island is introduced 
ClS db perturbation. 

The parameters used for these simulations are similar (but not equivalent) to those of the 



cyclone base case 42|. However, to simplify the physics we have set the temperature and 



density gradients in the background distribution {R/Lt = R/L^ = 0)to zero, otherwise 

• Inverse aspect ratio e = 0.19 

• Electron to ion temperature ratio T^/Ti = 1 

• Safety factor q = 1.5 and magnetic shear s = 0.16. 

• 2 toroidal modes, 167 radial modes. Results are presented with, k^pi = 0.025, k^pi 
0.05, fcjpi = 0.1 
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FIG. 6. (Color online) Normalized electrostatic potential {(pN = ecp/Tp^, W/pi = 24) in the plane 
perpendicular to the magnetic field (outboard mid-plane). Black lines represent the perturbed 
flux surfaces calculated from the total parallel vector potential. The top panel shows a vortex 
with a positive sign, while the lower panel shows that the vortex has flipped sign at a point later 
in the simulation. Data from Gyro-kinetic simulation without turbulence and initialised with a 
sin^ density perturbation which produces a electrostatic vortex structure similar to those seen in 
turbulent Gyro- kinetic simulations (See Fig. O) 



For an idea of the typical size of an island, consider a m = 3, n = 2 island which is 
resonant at g = 1.5. This choice effectively determines p^, = 4 ■ 10~^ in the ki^pi = 0.05 case 
and = 2 ■ 10~^ in the krPi = 0.025 case, values that correspond to a medium-size tokamak 



such as ASDEX Upgrade [43]. 



Run in this way, the code encapsulates all the necessary physics of the Geodesic acoustic 
mode in the MHD limit. A density perturbation is initialised which is resonant with the 
magnetic island, p = po sin ^ and then allowed to freely evolve. An oscillatory vortex mode 
is excited, as seen in Fig |6l The initial perturbation, however, is not an exact eigenfunction 
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of the system and thus we get some extra radial oscillations (See Fig. [6]) that damp away on 
a long time scale. 

Fig. [1] shows the time trace of the amplitude of the first non-zero electrostatic potential 
mode which has the same poloidal mode number as the magnetic island. The trace compares 
the potential from a gyrokinetic simulation (dashed line) without turbulence and one with 
(black line). 

From these traces we see that in both cases there is indeed two frequencies present. The 
first, faster component is the standard GAM oscillation produced by the Geodesic curvature 

38|. The second, 



of the circular flux surfaces which are also damped by kinetic effects | 
significantly longer period and higher amplitude oscillation, is the Geodesic mode around the 
closed flux surfaces within the magnetic island separatrix. This is slowly damped compared 
with the usual GAM case. There is indeed some disparity in the frequencies between the 
linear and nonlinear turbulent simulation traces, however the physics is significantly different 
between the two, with turbulence and the presence of equilibrium temperature and density 
gradients in the non-linear case, which could have a significant impact on the frequency of 
the vortex mode. 

Plotted in Fig. |5] is the squared oscillation frequency against the normalised squared 
toroidal wave- vector associated with the magnetic island, kg = m/r. We see that there 
is an exactly linear relation between these parameters. In the lower panel is plotted the 
squared frequency as a function of the flux surface label within the magnetic island as 
calculated from the eigenvalue solver with the same parameters as used in the gyro-kinetic 
simulations. Agreement between the simulations and the eigenvalue analysis is very good, 
with the frequency values matching near to the separatrix {Q = 0.95). An analysis with a 
kinetic or higher order theory which takes into account variation of the mode across the flux 
surfaces would give a more accurate calculation of the oscillation frequency, but it beyond 
the scope of this paper. 

It is observed that the mode from our eigenvalue analysis, with the closest matching 
frequency to that seen in gyro-kinetic turbulence simulations is the one which corresponds 
to an up-down density asymmetry within the magnetic island, with an eigenfunction corre- 
sponding to the full line in Fig. |3] and its symmetric reflection. This is depicted in cartoon 
form in the left panel of Fig. |H Solutions of this form have a GAM component, and therefore 
introduce a compression. The compression produces an electric field that is perpendicular 
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to the flux surfaces, and is tfierefore essential in producing the vortex structures as seen in 
Fig.! 

Indeed, solutions exist which have good agreement with between their eigenvalue and the 
frequency observed in gyrokinetic simulations, however, these have no GAM component and 
are therefore purely sound waves. Solutions of this form are unable to produce the vortex 
structures observed as plotted in Fig. [6] (and its electrostatic potential time trace in Fig. [1]). 
One example of this is the eigenfunction represented by the dot-dashed line in Fig. |3] whose 
eigenvalue is plotted (dotted line) in the lower panel of Fig. |5l 

In turbulence simulations it was observed that spreading occurred and turbulent struc- 
tures entered the island from the upper x-point and spread down into the island, giving an 
up down asymmetry (See Fig. 5 in Ref. IJ]). It is this mechanism which is a candidate to 
excite the oscillatory mode seen that has a value close to the calculated frequency near the 
separatrix and not a higher frequency as would be expected further toward the 0-point. 



VI. CONCLUSIONS 

Motivated by the observation of oscillatory potential vortex structures seen in simula- 
tions of turbulence around magnetic islands, we have performed an analysis of the Geodesic 
acoustic mode around a magnetic island. It is seen that long time-scale oscillatory solutions 
are generated with the same properties as the Geodesic acoustic mode, namely plasma com- 
pressibility producing an electric field perpendicular to the perturbed flux surfaces, which 
produce meso-scale vortex structures. 

Also performed were Gyro-kinetic simulations where the turbulence was neglected and 
the density initialised to be present with the magnetic island, which generates an oscillatory 
potential structure, which is the same as that seen in turbulence simulations. The scaling of 
the frequency of this oscillation agrees with our eigenvalue analysis, leading us to conclude 
that the oscillatory structures are indeed the Geodesic Acoustic Mode around the closed 
flux surfaces within the magnetic island. 

These oscillatory vortex structure can induce E x B flow around large magnetic islands 
which have a profound effect heat transport in the vicinity of the magnetic island and can 
also have a significant regulatory effect on the turbulence in this region. 
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VII. APPENDIX A - TREATMENT OF EXB TERM 

Firstly we consider the denominator of the compressive term, J JdS. If we utilise: 



VQ = ^ , + -smCe-- sm 34 
r R 



We utihse the magnetic field approximation B'^ — Bq/{1 + ecosOy and Taylor expand, 
to give: 

|2 



f^-^JdS^ (35) 



+ / JdS ^ ^/^ 4f^ + ^)(l + 2ecosg)sin^e 



We first perform the 9 integral, which removes the terms of order e, we obtain. 



27r / /"(fi + cos^) 2/"^ 



2 



16 / \^ , "^ d^ + I — + 

rfe). (36) 



9 \ /" • 2 ^ 
n \ I sm 4 



We assume a small aspect ratio (r < i?) and also small island width in relation to the 
minor radius (w < r) and as such the integral reduces to: 

Finally we treat the term / p BxvavB^ j^g^ utihsing B = + BeO + ^ sin Or and also: 



VE^ = 2BVB 



^ f cos - sin dd) (38) 



i?(l + ecose)3 
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Taking the cross product with the magnetic field vector: 



B X VB^ = BtC X VB^ + BeO x VB^ + VC x VV^ x VB' 

2R3 

° ' Bt cos 99 - Bg cos e( 



R{1 + e cos 



ihTTh 

+ Btsm9r - -—sm9sm^C) (39) 



Rr 

Taking the inner product with Vfi, this gives four terms: 

B X VB^ • = f- ^^^5. sin ^ 

R{l + ecos9y\ 

Btm . . 

cos y sm 4 cos 9 sm 4 



sin^^sin^) (40) 



— sir 

rR 

We make the approximations, neglecting the effect of the island on the field strength: 



Bt = fio/(l + ecos^) 

5, = ^/(1 + ecos^) (41) 



we finally obtain: 

BxVB^-Vn -2 /4(fi + cosO . 



B^ BoR V wy/2 



sin 9 



7TX €Tl 

H cos 6* sin 4' H — =: cos 9 sin ^ 

r git 

'"'''^sin2esin^(l + ecos^)) (42) 



rR 

The last term in Eq. fH2|) can be neglected as it is comparably small with respect to the 
other terms. 
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